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Experimental study of forces between quasi-two-dimensional emul- 
sion droplets near jamming 

Kenneth W. Desmond'', Pearl J. Young, Dandan Chen^, Eric R. Weeks* 
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We experimentally study the jamming of quasi-two-dimensional emulsions. Our experiments consist of oil-in- water emulsion 
droplets confined between two parallel plates. From the droplet outlines, we can determine the forces between every droplet pair 
to within 8% over a wide range of area fractions (j). We study three bidisperse samples that jam at area fractions 0c ^ 0.86. Our 
data show that for > the contact numbers and pressure have power-law dependence on — 0^ in agreement with the critical 
scaling found in numerical simulations. Furthermore, we see a link between the interparticle force law and the exponent for the 
pressure scaling, supporting prior computational observations. We also observe linear-like force chains (chains of large inter- 
droplet forces) that extend over 10 particle lengths, and examine the origin of their linearity. We find that the relative orientation 
of large force segments are random and that the tendency for force chains to be linear is not due to correlations in the direction of 
neighboring large forces, but instead occurs because the directions are biased towards being linear to balance the forces on each 
droplet. 



1 Introduction 

A liquid to amorphous- solid transition, also known as a jam- 
ming transition, occurs in a wide variety of soft materials such 
as colloids, emulsions, foams, and sand. In general the jam- 
ming transition is due to an increase in the particle concen- 
tration 0; the particles become sufficiently crowded so that 
microscopic rearrangements are unable to occur when exter- 
nal stresses are applied^ ^. At a critical 0^ the system jams 
into a rigid structure, and many of the material properties are 
known'SSto scale with a power-law dependence on {(j) — (j)c). 
' While these soft materials have obvious differences, it has 
been postulated that there are universal features of the jam- 
ming transition that all these materials share in common such 
as critical scaling and the emergence of force chains. 

In all systems above the jamming point, particles press into 
one another and deform. As the density increases, new con- 
tacts form and particles deform more, increasing the pressure. 
Interesting, both the average number of contacts z and the 
pressure P show critical-like scaling relative to the jamming 
point. In experiments and simulations, both 2D and 3D, the 
average number of contacts scales as z — Zc = A(0 — 0c)^S 
where Zc and A depend on the dimension and ft = 1 /2 re- 
gardless of dimensionEHHI. Simulations found P ^ (0 — 0^)^^, 
where fip depends on the details on the interparticle force 
If this pressure scaling and connection between be- 
tween and the interparticle force law extends to experi- 
ments, then this would demonstrate a direct link between the 
interaction of the constituent particles and the bulk properties 
of the sample, as the bulk modulus can be found from /^(0). 
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Another observed feature of jammed systems is the spa- 
tial heterogeneity of the particle-particle contact forces. In 
experiments and simulations, both 2D and 3D, the shape of 
the probability distribution of forces is broad with an expo- 
nential like tail^ ^ lo 21 ^j^^ largest forces tend to form chain 
structures that bear the majority of the load'^^tHl, These force 
chains are responsible for providing rigidity of jammed ma- 
terials to external st resses and are related to many other bulk 
propertiesEGHSEl. In prior experiments on 3D emulsions, 
the structure of the force chains was studied directly, where 
force chains extended over 10 particle diameters with an per- 
sistence length of 3 - 4 particle diametersEMEni^ 

There have been theoretical attempts to understand force 
chains, such as the q-model of Coppersmith et al^, directed- 
force chain networks of Socolar's group and simula- 
tionslSHlIl. Others took an ensemble approach to describe 
force chains, with different choices for ensembles^^ ^^ 
While some of these models successfully predict certain prop- 
erties of the force network, they can not explain the physical 
origins of force chains. To explain the structure of the fo rce 
chains observed in 3D emulsion studies, Brujic et al. Mand 
Zhou et al. proposed an accurate model that provides a 
physical description for the origin of force chains. This model 
has two simple assumptions: first, the forces on a droplet must 
balance, and second, forces between neighboring droplets are 
uncorrelated. This model has not been applied to 2D systems. 

In this paper, we introduce a new experimental system to 
study the universal nature of the jamming transition. Our sys- 
tem consist of quasi-2D soft deformable droplets with no static 
friction forces. In the appendix, we describe our method to de- 
termine the forces between droplets in contact to within 8%, 
significantly better than prior studies of foams and compa- 
rable to photoelastic disks Using our experimental model 
system, we find power-law scaling for the coordination num- 
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Fig. 1 (Color online) (a) A schematic of our co-flow apparatus. Oil 
is pumped at a constant rate through a micropipet centered within a 
capillary tube of larger diameter. Around the inner micropipet, a 5 
g/mL water Fairy soap mixture is pumped through the capillary 
tube, and as oil leaves the micropipet it forms spherical droplets that 
repeatedly break off with the same diameter, (b) A schematic of our 
sample chamber where emulsion droplets are confined to a 2D plane 
by two microscope slides separated by either a ~ 100 ^m spacer 
(transparency film) or ~ 180 ^m spacer (glass coverslip). 



ber and pressure (Sec. 4.2 ), we observe a relationship between 
the interparticle force law and j3p (Sec. 4.2 ), and see a distribu- 



tion of contact forces similar to prior work (Sec. 4.3 ). Further, 
we confirm the assumptions of the Brujic-Zhou model apply 
to our data and that the model well-describes our 2D data 



(Sec. |4.4[ ). This work provides an in depth study comparing 
data from our experimental model system to other numerical 
simulations, theory, and experimental systems, thus furthering 
our understanding of the jamming transition and supporting 
the applicability of ideas of jamming to a new system. 

2 Experimental method 

We produce emulsions using a standard co-flow micro-fluidic 
technique'^, see Fig.[T|a). The inner tube diameter is ~ 35 /im 
and the outer tube diameter is ~ 500 /im. The continuous 
phase is a mixture of water and the commercial soap "Fairy", 
and flows through the outer tube at a rate of ^ 1 mL/min. The 
droplets are mineral oil, which flows through the inner tube at 
a rate of ~ 0.5 /iL/hr. Slight variations of these parameters let 
us produce monodisperse droplets with radii in the range of 
80-170 /im; any given batch of droplets has a polydispersity 
of less than 4%. Mixing together two monodisperse batches 
lets us produce bidisperse samples with whatever size ratio 
and number ratio desired. 

Our sample chamber is designed to create a system of quasi- 
2D emulsion droplets, analogous to 2D granular systems of 
photoelastic disks but without static friction. The chamber 
consists of two microscope glass slides of dimensions 25 mm 
X 75 mm (Coming) separated by a ~ 100 {im spacer (trans- 
parency film) or a ~ 180 {im spacer (Corning No. 1 glass 
coverslip) glued along the two longer edges; see Fig. [ijb). 
The sample chamber thickness is tuned so that the droplets 
are deformed into pancake shapes, with aspect ratio (diame- 
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ter/height) ranging from 1.6 to 3.0; see Fig. 12 



Fig. 2 (Color online). Illustration of our experiment. Oil droplets 
rise to the top of the sample chamber due to buoyancy. At the bottom 
of the droplet "pile," droplets barely touch and are not deformed. At 
the top, droplets are compressed due to the buoyant weight of the 
droplets below them. This lets us study the sample from the 
jamming area fraction on up, and also provides a means to calibrate 
the forces as described in the Appendix. The scale bar is 200 jUm, 
and the images have area fractions = 0.88, 0.92, and 0.96. 



After the sample chambers are filled, they are placed on 
a microscope for imaging with either a 1.6x or 5x objec- 
tive lens. The droplets are allowed to equilibrate their po- 
sitions; we only consider static samples. Our camera takes 
2, 200 X 1, 800 pixel^ images. We overlap images from differ- 
ent areas to construct a single large field of view image on the 
order of 10,000 x 50,000 pixel^ containing between 1,000 to 
5,000 droplets depending on the droplet sizes. We image every 
droplet (wall-to-wall) and only analyze droplets more than 4 
diameters away from the nearest wall to avoid wall effects . 

3 Empirical Force Law 

We wish to use the droplet images (such as those shown in 
Fig. [2]) to determine the forces droplets exert on each other. 
An isolated droplet is circular with 2D radius ro due to sur- 
face tension. Droplets feeling forces from other droplets are 
deformed. Our goal is not to know the exact form of the force 
law governing inter-droplet forces. Rather, we need to know 
these forces to within our experimental error. The details of 
our approach are given in the Appendix; we briefly summa- 
rize our method here. 

By tilting the sample as shown in Fig. [2j we exploit the 
known buoyant forces on the droplets, which are O(10~^ /iN) 
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Fig. 3 Close up view of regions within our three bidisperse samples. 
The scale bar in each image is 250 jUm. (a) is an image for our 
sample with size ratio 1.25 at = 0.89. (b) is an image for our 
sample with size ratio 1.42 at = 0.90. (c) is an image for our 
sample with size ratio 1.52 at = 0.90. 



per droplet. Droplets further up the incline feel larger forces 
and deform more because they must support the buoyant 
weight of the droplets below. This can be seen in Fig.[2j where 
droplets further up the incline are more compressed. 

The key parameters to determine the forces are the contact 
length A/ of the interface shared between the two droplets, and 
the modified radius of curvature r of each droplet deduced 
from the perimeters of the droplets where they are not con- 
tacting another droplet. Laplace's Law relates the interface 
curvature to the pressure difference between the interior and 
exterior of a droplet; in 3D this relation is AP = y/r with y 
being the surface tension. From this, it is clear that deform- 
ing a droplet (increasing its internal pressure) results in r < tq. 
Likewise, it is clear from Fig. |2] that A/ is larger for droplets 
feeling larger forces (droplets farther up the incline). 

To find a force law from our data, we use numeric methods 
to find a force law /(A/,r) that best satisfies what we know 
about the data: (1) the net force on any droplet is zero (as the 
droplets are motionless), (2) the components of the forces in 
the y direction balance the buoyant weight of each droplet, and 
(3) the forces acting between two touching droplets are equal 
and opposite. After some work, we find for two droplets / and 
7 in contact that the best functional form is 

F = ai (R^Llnj) + a2{RoL/rijf, (1) 

where l/vij = l/r^ + l/r^, Rq is the mean 3D droplet radius 
prior to putting the sample into the quasi-2D chamber, and ai 
and a2 are constants that depend on the oil- water- surfactant 
details. For our samples, ai = 4.25 /iN/mm and a2 = 4.12 
/iN/mm^. The forces we observe in the experiment range 
from - 0.5/iN, several hundred times larger than the buoyant 
weight of an individual droplet, which makes sense given that 
the inter-droplet forces at high (j) are due to the accumulated 
pressure from the buoyant forces of the many droplets below. 
Intriguingly, we find that the force law does not depend on the 
gap thickness of the sample chamber. We stress that Eqn. [T] 
is an empirical deduction and is only approximate. The data 
presented in the next section are from images taken with the 
L6x objective lens and the forces obtained have a random un- 



G Hr r^'^g^ poly(^^g) r^"^'"'^ pol^^ 

1.25 0.684 126 Aim 3.4% 102 ^um 3.1% 

1.42 0.849 130 Aim 3.0% 105 Aim 3.4% 

1.52 0.806 137 3.3% 90.1 Aim 3.1% 



Table 1 This table provides parameters characterizing the droplets 
in our 3 bidisperse samples. The first column o = r^^^"^ /tq^^^^"^ is 
the size ratio and the second column is the number ratio 
rir = A^(^^g) /A^(small) polydispersity in droplet sizes for the big 
droplets is indicated by poly*^^^^) and for the small droplets by 
poly(s"^^ii). jj^g polydispersity is defined as the standard deviation in 
droplet sizes normalized by the mean size. Our uncertainties are ± 2 
Aim in tq, ± 0.1% in a, and < 0.1% in polydispersity. Since we 
image every single droplet there is no measurement error in rir. 

certainty of 16%. For a more detailed discussion on obtaining 
this force law and the uncertainty, see the Appendix. 

4 Jamming of Binary Packings 

We will investigate the jamming transition and force networks 
in disordered bidisperse packings using the setup shown in 
Fig.|2] By allowing the sample to equilibrate in this chamber, 
we set up a gentle gradient in area fraction ranging from just 
barely jammed at the bottom of the incline to well-jammed 
near the top. This lets us study the properties of the pack- 
ing over this full range of area fractions. For the analysis 
in this section, we only consider droplets between 5 and 50 
droplet diameters up the incline, unless otherwise specified. 
This region corresponds to an area fraction (j) ranging from 
0.89 — 0.96. We study three different samples with different 
size ratios, as shown in Fig. [3] Details of each sample are 
given in Table [T] 

4.1 Identifying jamming area fraction (j)c 

We start by identifying the jamming area fraction (j)c for each 
data set. To determine the jamming point, we calculate the 
area fraction with distance up the incline directly from our 
images of the droplets in three steps. First, using the cen- 
ters and radii of the droplets we compute the radical Voronoi 
cellsM 

for each droplet. The radical Voronoi tessellation 
divides space into polygons, one per droplet, taking into ac- 
count each droplet 's size so that each droplet is fully con- 
tained within its own polygon. Second, we determine the 
area Ay of each Voronoi cell and the area of each droplet. 
Third, the area fraction (j){y) at a position y is computed as 
0(j) = Zy^Aj y^/Zy^Ay y^, whcrc k indcxcs all droplets with a cen- 
ter of mass within y — Ay/2 and y-\-Ay/2. For this step, all 
the droplets are examined, down to nearly y = 0, except for 
the droplets at j ^ where the Voronoi cell is poorly defined. 
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We choose Ay = 6(ro) where (ro) is the mean droplet radius 
(which depends on the size ratio and number ratio, see Table 
[T]for details). This value of Ay gives roughly 150 droplets per 
y sampled. Within this window of Ay, = {d(l)/dy)Ay ^ 
0.007. From 0(_y) we can obtain the jamming point 0^ by ex- 
trapolating the value of to j = 0, where }7 = is defined 
as the bottom of the droplet pile. We can treat the j = 
point in our data as the jamming point since the forces be- 
tween droplets at y = are nearly zero. For the three data 
sets, the extrapolation is done by fitting (j){y) to a power law 
[(l)(y) = 0^ + giving (l)c =0.855 ± 0.005, 0.861 ± 0.005, 
and 0.858 ± 0.008 for the data with size ratio a = 1.25, 1.42, 
and 1.52, respectively. We chose to use 0(j) = (l)c-\-ay^ since 
(j) — (l)c vs y appears fairly linear on a log-log plot. In simu- 
lations on frictionless disks and experiments on 2D f oams it 
has been reported that 0c ^ 0-84 for bidisperse systemsl^CIEl^ 
which is a little lower than the values we found. 

Our measured area fraction depends on where we define the 
outer perimeter of a droplet. As seen in Fig. [3] the droplets 
have thick black outlines. We look at the outer edge of each 
outline, and define the perimeter as the pixel location where 
the intensity is halfway between the white color outside the 
droplet, and the black color in the darkest part of the outline. 
The transition from black to white occurs over a distance of 2- 
3 pixels, and so we judge that we have a systematic uncertainty 
in the area fraction of roughly 1% due to the determination of 
the perimeter position. Since this is systematic, the distance 
to the jamming point ((j) — (j)c) is insensitive to this error and 
therefore in most of our results we focus on — 0^. 

4.2 Critical Scaling 

The first critical scaling we investigate is the coordination 
number, the mean number of contacts each droplet has. Prior 
numerical studies of jamming in frictionless systems found 
that the coordination number z obeys a power law scaling of 
the form z — Zc = A(0 — 0c)^S where A ~ 3.5, Zc = 4, and 

= l/2'3'3. It has been observed that A has a slight de- 
pendence on the force law and polydispersity, but Zc = 4 and 
P = 1 /2 are independent of the force law and polydispersity. 
Katgert et al ^ ^ found for a 2D bidisperse foam with size ratio 
1.5 a critical scaling with A = 4.02 ±0.02 and ft = 0.50 ±0.02 
while fixing = 4. The critical point Zc has been interpreted 
as the isostatic point Ziso (minimum number of contacts nec- 
essary for a mechanically stable packing). For 2D, Ziso = 4, in 
agreement with Zc found in prior work. 

To compare experimental data and simulation data, the ex- 
perimental area fraction needs to be converted into a theoreti- 
cal area fraction ^ ^ . This is because the simulated particles are 
allowed to overlap (thus diminishing the total area they take 
up at large (j)) while our experimental droplets always occupy 
the same total area. We convert our experimental (j) values to 
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Fig. 4 (Color online) (a) Scatter plot of coordination number against 
^theory — 0c- All data were fitted together to z — = A(0 — 
where the fit is shown as the black dashed line with fit parameters 
Zc = 4.2, A = 3.2, and = 0.4. Fitting the different data sets 
separately gives slightly different fit values, listed in Table|2] (b) A 
scatter plot between pressure and ^theory — 0c- The pressure has been 
scaled by c = 1, \/T0, and 10 for the G = 1.25, 1.42, and 1.52 data 
respectively. Each data set is fitted to ^ = A(0theory — 0c)^^» shown 
as the black dashed lines. The fit values are given in Table |2] 



^theory using the method of Katgert et al^^. From our data we 
determine z and ^theory at various points along the incline. The 
results are plotted relative to the jamming point in Fig. ^sl), 
and show power-law scaling. Fitting the each data set to the 
theoretical scaling law, z — Zc = A{(j) — (j)c)^^,wQ obtain values 
for A, Zc, and jS^ which are reported in Table [2] Our values of 
A ^ 3.2 are close to A ^ 3.5 found in a numerical study by 
O'Hern et al ^ for particles with size ratio 1.4. The fitted val- 
ues for Zc are w ithin the uncertainty of the previously found 
value of 4'^^'^^'. However, our droplets have a slight attraction 
which may result in a slightly tighter packing of droplets at 0^ 
with a coordination number > 4. Given our uncertainties of 
Zc, our data are consistent with both Zc = 4 and Zc > 4. Finally, 
for each packing, the exponent ^ 0.4 agrees with the prior 
findings (j3 = 0.5) to within our uncertainty, although we have 
a fairly large uncertainty in our exponents. Interestingly, in 2D 
photoelastic disk experiments, they found z — Zc = {<\> — <\>c)^^ 
with = 0.53 ± 0.03 without needing to convert their exper- 
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2.4 ±0.1 1.19 ±0.02 




1.52 


2.0 ±0.1 1.15 ±0.03 





Table 2 The fitting parameters for the power law fits to the data for 
each size ratio o. Note that simulations found j3p = jSyt^ see text for 
a discussion. The uncertainties in the fit values are obtained by 
computing the standard error in each fitting parameter. 



imental to ^theory but A ~ 25 for that study which is con- 
siderably different from our results. In their work, they were 
limited to area fractions close to 0^ due to the difficulty of 
compressing their particles to high area fractions, while our 
data (and those of Ref . ^ ^ ) extend over a larger range of . 

The second critical scaling we investigate is the dependence 
of pressure with distance to the jamming point. Simula- 
tions of 2D particles found = A(0theory — ^c)^^ ^ where A 
and j3p depend on the form of the force law. In the numeri- 
cal study by O'Hern et alW, they used frictionless disks that 
interacted via the force law fij = Fo{5rij/dij)^f , where Fq is 
a scale, 5^^- is the distance between two particles in contact, 
and dij is the sum of the radii of the particles in contact. They 
found that jip = jif. It is certainly possible that for other force 
laws, the scaling of pressure with {(j) — (j)c) could be differ- 
ent. In particular, in our experiment, the force between two 
droplets is not a unique function of Srij but rather depends 
on the droplet perimeters which are influenced by all of their 
neighbors. In 2D photoelastic disk experiments jip was found 
to be 1.1 No prior experimental 2D studies have examined 
the scaling of ^ for systems without static friction. 

For our experiment, we compute the local pressure of our 
sample by first locating a set of droplets k within a window 
y — Ay/2 and y-\-Ay/2. For these k droplets the pressure is 
^ = Y.ilLj>i^ij^ij l^k^v^ where / and j index all contacts on 
the k droplets and ^k^k.v is the sum of the Voronoi areas of all 
k droplets^^. In this formula, Fij and Vij are both taken to be 
positive scalars. Here we use Ay = 5ro. In Fig. |4|b) we plot 
the pressure for all three packings against (^theory — 0c- These 
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Fig. 5 (Color online) The average force between droplets in contact 
plotted against the amount of compression between the droplets. 
The average force has been scaled by a prefactor of c = 1 , >/T0, and 
10 for the a =1.25, 1.42, and 1.52 data respectively. Each data is 
fitted to (fij) = Fo{Srij/dij)Pf and the fits are shown as the black 
dashed lines. The fit values are given in Table [2] Note that this data 
is an effective force law, not the true force law: for a given 5rij/dij, 
different droplet pairs may experience different contact forces. To 
illustrate this, we have added error bars to the plot, where the error 
bars represent one standard deviation in the spread of measured 
contact forces at each Srij/dij. 



results show power-law scaling. The dashed lines are the fit to 
^ = A ((^theory — ^c)^^ with the fit valucs shown in Table [5 In 
particular, we find j3p values between 1.26 - 1.41, larger than 
j3p = 1.1 found for photoelastic disks 

To compare with the simulations of O'Hern et alW, we wish 
to approximate how forces between our droplets depend on 
their separations 5^^. For each observed we find the true 
force fij from our force law. We average all of the observa- 
tions over small windows in to find an effective average 
force law as a function of 5^^, plotted in Fig. [5] The error 
bars emphasize that Fig.|5]is only an average trend rather than 
the true force law. Intriguingly, the averaged data follow a 
power law: we fit each data to {fij) = Fo{5rij/dij)^f to obtain 
the power law exponent The fits are shown as the black 
dashed lines in the figure, with fit values listed in Table [2] 

Our fits give j3/ < j3p in contrast to the results of O'Hern 
et ah ^ where j3p = j3/. This equality was found for systems 
close to the jamming area fraction. The exponent for the pres- 
sure, relates to how droplets are compacted with increas- 
ing (j) — (j)c^. Close to when (j) is slightly increased, droplets 
can avoid significant compression by rearranging and form- 
ing more contacts, however, at larger 0, droplets can not form 
many new contacts and must instead undergo larger compres- 
sion. Therefore, at larger area fractions, the pressure increases 
more rapidly with (j) — (j)c than it does near the jamming point. 
This argument predicts Pp > Pf, in agreement with our data 
which extends far from 0^- While the uncertainty in each force 
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measurement is 16%, this uncertainty is unlikely to signifi- 
cantly affect the pressure results, as the data of Figs. [4] and [5] 
are averages over many forces. 

4.3 Force Distribution 

We now consider the distribution of contact forces for each 
packing at different area fractions. Like before, we sample the 
contacts forces at various points up the incline using a window 
ofy — and However, we need many contacts 

to obtain a good distribution of contact forces, and therefore, 
we use = 30ro. (Over this range of A_y and for droplets at 
least 10 diameters up the incline, = {d(^/dy)^y = 0.025.) 
This window size gives roughly 2500 contacts for each y sam- 
pled. In Fig. |6ja), the lines show the distribution of contact 
forces normalized by the mean contact force at locations with 
— 0c as indicated; all the data are for the a = 1.25 pack- 
ing. In all our data most forces are near or less than the mean 
force (/) and that the maximum force is about 3(/), with a 
somewhat exponential tail. The shape and magnitude of all 
the curves are roughly the same. All curves show a dip at 
small forces. The symbols in Fig. |6ja) show the distribution 
of normal forces from Majmudar et ah an experiment using 
frictional 2D photoelastic disks. In their experiment, the parti- 
cles were isotropically compressed to an area fraction ~ 0.016 
above the critical area fraction. Our results look essentially 
the same as theirs, despite the differences in experimental sys- 
tems. These distributions are also similar to simulations and 
experime nts on 2D and 3D frictionless systems on foams and 
emulsionsl2inC2IIIl, Xhis suggest that the distribution of forces 
is a fairly universal property of all jammed systems. The shape 
of our measured distribution is also in general agreement with 
the predictions of the force network ensemble and is also in 
agreement with the earlier q-model^^, although our data are 
not clear enough to distinguish between these two models. 

As observed by Katgert et al.^^, our force distributions 
slightly narrow with increasing area fraction in Fig. [6ja). To 
quantify this, we plot the standard deviation of each force dis- 
tribution in Fig.|6jb). For the cr = 1 .25 and cr = 1 .42 packings, 
the width of the force distributions decreases quickly above 
the jamming point and then saturates to a constant width of ~ 
0.5. The O" = 1.52 packing has a broader distribution of forces 
at larger — 0c compared to the other two packings, and the 
decrease in the width as — 0^ increases is more subtle. Over- 
all, our results are qualitatively in agreement with ^ ^ , although 
they did not compute the standard deviations. It is possible 
that measurement errors in our forces have a dependence on 0, 
as discussed in the Appendix, and that this could affect the de- 
pendence of the standard deviation on (j) . However, the errors 
are no worse than 16%, much less than the width of the distri- 
butions for any 0, and so a possible -dependence of the errors 
cannot account for the decrease in width seen in Fig.[6jb). 
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Fig. 6 (Color online) (a) Distribution of contact forces relative to 
the mean contact force at different (j) for the a = 1.25 packing. The 
mean force is (/) = 0.011 jdN, 0.045 fdN, and 0.13 j^N for the 
(^>-(^>c = 0.024,0.062, and 0.106 data respectively The solid 
triangles are data from a 2D photoelastic disk experiment at 
— 0c ^ 0.016 from Majmudar et ah (b) The standard deviation 
of P{f / (/)) for each packing at different 0. The standard deviation 
of the Majmudar et al data is 0.52. 



4.4 Force Chains 

In this section we consider various statistical measurements 
on the randomness of the force chain network, and then we 
compare predictions of the Brujic-Zhou model to our 

data. When analyzing data in this section we consider all 
droplets and contacts between 40 < j/2ro < 80, and over this 
range increases from 0.93 to 0.96. We find that all of the 
properties discussed below do not depend on ^ at larger area 
fractions, and so considering this larger range of gives us 
better statistics. 

To start, we define a force segment to belong to a force chain 
if it is one of the two largest forces on both droplets joined by 
the force segment. Under this definition, each droplet can only 
have a maximum of two force segments that belong to a force 
chain, and therefore, our definition does not allow for force 
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Fig. 7 (Color online). This image shows only the forces belonging 
to a force chain within a region of the a = 1.25 sample. On average, 
the forces are larger further up the image because the sample is 
inclined, and this can be observed in the image by the increasing 
redness of the force segments at the top. 




Fig. 8 (Color online). Definitions of the angles 6i and 62 between 
joining force segment. In the sketch both 61 (clockwise to extended 
line) and O2 (counter clockwise to extended line) are positive. If 
there is a correlation in orientation that tends to make force chains 
linear, then the correlation between 61 and 62 is positive. 

chain branching or merging, which will simplify the analysis 
below. Note that even droplets with small forces can partici- 
pate in force chains, given that the definition only requires the 
force segments to be large for the given droplet and its con- 
tacting neighbors. Figure [7] shows force chains using our def- 
inition. These chains are fairly linear and vary in length from 
1 to more than 10 force segments. Our samples are formed by 
slow compaction due to buoyant forces, which m ay introduce 
an anisotropy in the force chain networlJSIllEII. Indeed, we 
find a slight bias for force chain segments to be orientated up 
the incline. 

Given that force chains form linear like structures and that 
there is a slight anisotropy for force chains to align along the 
incline we may expect there to be correlations in the orienta- 
tion of neighboring force chain segments. To quantify such a 
correlation we define two relative angles 61 and 62 between 
joining force segments, where the definition of 61 and 62 are 
shown in Fig. [8] We compute the Pearson correlation coeffi- 
cient C = cov(0i ^62)1^1^ where cov(0i , O2) is the covariance 



of Q\ and 02 and Gq is the standard deviation of Q. We find 
that for all C is zero or nearly zero (at most C = 0.2), in- 
dicating no cor relation. This agrees with prior work on 3D 
emulsions Thus, the apparent linearity of force chains 
seen in some locations of Fig. |7]is not due to correlations in 
the relative direction of neighboring segments that would keep 
the chain straight. 

To further explore the tendency for force chains to be lin- 
ear, we consider the distribution of Oi , where we drop the sub- 
script 1 as we are only focusing on two force segments at a 
time rather than three. In Fig. [9ja) we plot the distribution 
in 6 for all three packings. The distribution shows that most 
force chain segments form at an angle \ 0\ < 60°. Thus, force 
segments tend to form a linear chain not because their orien- 
tations are correlated, but simply because it's more probable 
that they are oriented at small angles relative to each other. 
Using our P{6) data, we determine the persistence lengths / 
using the standard definition of persistence length for polymer 
chains. 

We find / = 4.4(ro), 4.8(ro), and 3.8(ro) for the a = 
1.25,1.42, and 1.52 data. These are the distances beyond 
which the force chain has "forgotten" its original direction. In 
analyzing the distributions similar to P{G) for 3D emulsions, 
Zhou et al found a persistence length slight larger around 
/-6-8(ro). 

To further consider the orientations of force segments in 
force chains, we consider a model proposed by Brujic et 
^/ li9|20| extended by Zhou et al}^^. The Brujic-Zhou 
model is a method for generating ensembles of local particle 
configurations (a central particle and contacting first neigh- 
bors) and the forces acting on a central particle by its first 
neighbors. Each local configuration is generated by randomly 
placing Zi contacting neighbors such that any two neighboring 
particles do not overlap. Next, the contact forces between the 
central particle and Zi — 2 neighboring particles are chosen at 
random from a distribution P{f), leaving two unknown con- 
tact forces. We choose P{f) to match our experimentally mea- 
sured distributions (see Fig.|6]). By invoking force balance, the 
two remaining contact forces are found algebraically. Once 
a sufficient number of local configurations are generated, the 
distribution of force chain orientations can be studied. The ba- 
sic assumptions of this model are force balance, randomness 
in the magnitude of forces, and randomness in the orientation 
of forces. For our data the first assumption applies because 
the system is in mechanical equilibrium and above we have 
shown that the other two assumption reasonably apply. 

One issue in using the Brujic-Zhou model to predict is 
that the model only gives the forces between a central droplet 
and its first neighbors. To define a force chain segment we 
also need to know all the forces acting on each first neigh- 
bor as well. We therefore extend their model by generating 
additional forces on the neighboring droplets in exactly the 
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Fig. 9 (Color online) (a) Distribution of 6 for each packing, where 
both 6i and 62 are treated as a single variable 6. The red solid line 
is the distribution for the a = 1.25 packing, the green dashed line is 
the distribution for the a = 1.42 packing, and the blue dashed-dot 
line is the distribution for the CJ = 1.52 packing, (b-d) Comparisons 
between the experimental distributions and the predictions of the 
Brujic-Zhou model, for size ratios (b) a = 1.25, (c) a = 1.42, and 
(d) CT = 1.52. 



same way (constrained by the forces already chosen for the 
central droplet). This lets us apply our force chain definition 
given above, which requires that force segments be among the 
largest two forces on both droplets the force acts between. 
We repeat this extended Brujic-Zhou algorithm many times 
to compile data from all cases where the algorithm gives an 
instance of two valid force segments so that we can determine 
61 . To make the inputs into the model as consistent as possible 
with our experimental data, instead of randomly generating 
local configurations, we randomly select local configurations 
from our experimental data. 

Figure [9jb)-(d) compares P{6) measured in our experi- 
ments (black solid curves) with the predictions of the model 
(red dot-dashed curves). The model is in good agreement with 
the experiment, with the exception of some discrepancies in 
the magnitudes of the peaks. The model captures significant 
features of the data: for instance, the peak around = 0° 
is much different between Fig. [9jb) and Fig. [9jd), and the 
model replicates this difference. We also note that if we loosen 
the definition of force chain segments to simply those forces 
that are the largest two forces acting on any droplet (indepen- 
dent of how large they are relative to forces on neighboring 
droplets), we find nearly identical distributions as the ones 
shown in Fig. [9] 

Our analysis suggests so far that the force chain network is 
random, without long-range correlations. It therefore seems 
plausible that the distribution of force chain lengths should 
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Fig. 10 (Color online) Distribution of the number of force segments 
making up distinct force chains. The data points are experimental 
values and the dashed lines are fits to the data of the form 
P{n) = (1 where p is found to be 0.722, 0.758, and 0.717 

for the (7 = 1.25, 1.42, and 1.52 packings, respectively. 



obey a random process. If there is a probability p for a force 
chain segment to be connected to a neighboring force chain 
segment, then the distribution of chain lengths should obey 
the scaling P{n) = {I— p)p^, where n is the number of force 
segments within a force chain. 



In Fig. [T0| we plot the dis- 
tribution of chain lengths for each packing. The data decay 
exponentially over 3 orders of magnitude. The data are fit by 
P{n) = {l—p)p^ with p ^ 0.73 (see caption for details), indi- 
cating that it is highly likely that for a force chain to propagate 
through the material. The fits are shown as the dashed lines 
and show good agreement with the data other than at ^ = 1 . 

In granular quasi- static intruder simulations with friction 
between particles by Peters et ah , using a more sophisticated 
definition of force chains, they also found an exponential dis- 
tribution of chain lengths From their reported data on P{n), 
we estimate a value of = 0.65. It appears that statistically a 
force chain can be thought of as a random process with proba- 
bility p for the force chain to propagate, independent of but 
perhaps depending on the sample details. 

5 Conclusions 

We have introduced a new experimental model system com- 
posed of quasi-2D emulsions droplets to study the jamming 
transition. Our droplets are circular in shape and deform 
when press into one another, and at the contacts between 
two droplets the forces are in-plane mimicking a true 2D sys- 
tem. We can accurately measure the forces between touching 
droplets to within 8%, where our method is not limited to our 
experiment, and could be extended to determine forces in 2D 
foams, 3D emulsions, and 3D foams. Our model system has 
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unique strengths; we can easily make samples with any distri- 
bution in particle sizes, emulsions are stable over many days, 
setup is cheap, our droplets have no static friction, and our 
method can be extended to cases of flow"^^. 

Using our model system we observed power-law scaling of 
the contact number and pressure with — 0c, similar to prior 
numerical models^ ^. Notably we find that all three fit pa- 
rameters for the contact number scaling are quite close to the 
values found in 2D simulations. We verify experimentally for 
the first time a link between the interparticle force law and 
the critical pressure exponent, illustrating a direct relationship 
between the bulk properties of an amorphous solid and the 
interaction between the constituent particles. The agreement 
of our results and the numerical models shows that the qual- 
itatively different particle interaction we have does not play 
a significant role in determining the geometric structure and 
bulk modulus. 

Our analysis of the inter-particle forces found a probabil- 
ity distribution of forces in good agreement with those found 
in prior experiments and simulations, strongly suggesting that 
the shape is universal. We further examined the spatial struc- 
ture of the large forces ("force chains"). The directions of 
neighboring force chain segments are uncorrelated although 
there is a tendency for two force chains to be in the same di- 
rection. This is a sensible result as this allows the large forces 
acting on a droplet to balance one another. The Brujic-Zhou 
model, which assumes random and uncorrelated force seg- 
ments, recovers our experimentally observed probability dis- 
tribution of angles between adjacent force segments. 

This work provides more evidence for the universality of 
various properties of the jamming transition, such as critical 
scaling, the shape of the force distribution, and the structure 
of the force network. 
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A Method For Determining Force Law 

In this section we describe in detail our method for determin- 
ing an empirical force law that relates the outline of droplets 
to the contact forces. For an overview of our method see 
Sec. |3] This section is organized in the following manner: 
first, we discuss the measurements from droplet outlines; sec- 
ond, we discuss the general form of possible force laws; third. 




Fig. 11 (Color online), (a) Experimental image of droplets with the 
contacts found using our algorithm indicated by green lines, (b) 
Experimental image of a rattler droplet being held in contact with 
two neighboring droplets due to a slight adhesion. In this image the 
droplet are motionless and the system is not inclined, (c, d) A close 
up view of a droplet. Each portion of the water-oil interface is fitted 
to an arc with constant radius of curvature. The fits are shown as the 
different colored arcs. 



we present the optimization problem; fourth, we deduce the 
best force law consistent with the data. 



A.l Measurable Variables 

In this subsection, we discuss the various quantities measur- 
able from droplet images, and their measurement errors. In 
the following subsections, these quantities will be used to de- 
termine the forces between droplet pairs. 

The larger the contact between two droplets, the more force 
they feel. This is quantified by the contact length between 
two droplets. We measure this by identifying the portion of 
each droplet' s perimeter that is shared between them, shown 



as the light green lines in Fig. 11a). Uj is calculated as the 
length of the line segment. Since we can only measure the two 
endpoints of each contact to 1 pixel accuracy, we have an un- 
certainty 51 of >/2 pixels. For our highest magnification lens 
(5x) this gives 5/ = 1.1 jim and for our lower magnification 
lens (1.6 x) this gives 51 = 1.96 /im. 

One expects that any two droplets with a nonzero con- 
tact length Qij > 0) would experience a repulsive force at the 
contact. However, we observe a slight attractive interaction 
between droplets as shown in Fig. [TTJb), where the central 
droplet is adhering to two neighboring droplets. In this par- 
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ticular configuration, the adhesion force is balanced by the re- 
pulsive force, and therefore, the net force at each contact on 
the central droplet is zero for some finite contact length /q. 

We determine /q by averaging together the contact lengths 
found between 25-50 droplets at the very bottom of the in- 
cline. The standard deviation of these contact lengths is about 
0.025 ro, which is one to two orders of magnitude smaller than 
measured contact lengths lij of deformed droplets in contact, 
and therefore we conclude that /q is reasonably well-defined 
by the mean value. In addition to the adhesion, some of the 
finite contact length is due to optical resolution limits result- 
ing in a systematic effect: for similar conditions, /q is shorter 
when using the 1.6 x lens (data in Table 3) as compared to the 
5x lens (data in Table 4). Accordingly, the force law will be 
taken to depend on A/ = / — /q, canceling the systematic ef- 
fect, and the force law will be required to obey /(A/ = 0) = 0. 
We could also require /(/) — /(/o), but this will not produce 
force laws of any better quality and often the force laws will 
be more complex in form. 

Near the jamming area fraction, our fractional uncertainty 
is typically 5l/Alij ^ 100$. In contrast, ten droplet diameters 
up the incline ((j) — (j)c ^ 0.01) typical values of Alij are much 
larger and the uncertainty drops to about 51 /Alij ^5%, and 
decreases further still for larger (j) . 

Next, we wish to know the curvature of the droplet outline. 
We measure n for a droplet by locating each portion of the 
perimeter belonging to a water-oil interface, that is, not touch- 
ing another droplet. We fit these portions of the interface to 
an arc of constant curvature as shown in Fig. [TTJc,d) to ob- 
tain a local radius of curvature for each portion. The radii for 
the different portions should be the same, but in practice they 
vary due to noise. We arithmetically average these local radii 
of curvature to obtain the mean curvature for the droplet. To 
determine the uncertainty, we create artificial images of cir- 
cular perimeters with known radii and noise commensurate to 
our experimental data, and fit these perimeters to find values 
of r. The distributions of r reproduce the experimentally ob- 
served variance in the individual radii of curvature, and let us 
deduce that the measurement error of is 3.5% for the 5x 
lens and 7% for the 1.6x lens. In contrast to Al which is eas- 
ier to measure at large 0, is harder to measure at large (j) as 
the curved portions of droplets are shorter. Closer to these 
uncertainties are smaller. 

The last measurable quantity to consider is R, the mean 3D 
curvature of a compressed emulsion droplet in our quasi-2D 
system. This curvature relates to the Laplace pressure and so 
may be relevant for the force law, although we will show be- 
low that it is not needed; nonetheless we discuss it for com- 
pleteness. For scenarios where droplets are asymmetrically 
deformed in 3D, the water-oil interface has two principle radii, 
the maximum radius of curvature Ri i and the minimum radius 
of curvature Ri 2- For droplets compressed in this manner, the 




Fig. 12 (Color online). An experimental image of a mineral oil 
droplet squeezed between two glass slides, where the gap thickness 
is 1 mm, 7?^ II = 0.88 mm and = 0.56 mm. The orange (light) 
dashed line is a fit to the perimeter to obtain 7?^ y and Rf j^. 



mean curvature I / Ri = 1/2(1/ Ri^i -\- l/Ri^i) is constant any- 
where on the surface. 

To measure Ri \ and Ri 2 experimentally we take side view 
images of isolated droplets in a sample chamber of gap thick- 
ness /z = 1 mm (see Fig. [12]). The width of the droplet cross- 
section is 2Ri \, and corresponds to the droplet radius that 
would be measured as the 2D radius in the normal top-down 
view of our experiments. The free surface of this compressed 
droplet is a surface of mean curvature R\ this is not a circular 
arc of constant radius as R\ varies with height. To obtain Ri \ 
and Ri 2 of the droplet, we fit the surface using the the method 
of Caboussat and Glowinski^^ (an algorithm to generate the 
surface of a droplet compressed between two boundaries). In 



Fig.[T2j we show the fit as the orange dashed line. Repeating 
this method for many droplets, we find 7?;^//z = 0.552±0.011 
for droplets in the size range we use. For simplicity, we simply 
use Ri^^ = 0.552/z for all r,. 

A.2 Mathematical Treatment of an Empirical Force Law 

Our goal is to find an empirical force law /(//j,/o,r/,rj) re- 
lating the contact force between two droplets / and j to the 
information about their outlines. A priori it is useful to con- 
sider what such a force law should look like. 

We first consider two cases where the force law is already 
known, the ideal 2D case and the ideal 3D case. By ideal, we 
mean that the contact angle between two droplets is zero, and 
where there are no adhesive forces. Generally these are not 
realistic assumptions, due to the interactio ns bet ween the sur- 
factant molecules at the contacting interfaceEEl, por the ideal 
cases, the force between two droplets in contact can be mod- 
eled using Princen's 2D model ^^^^ or Zhou's 3D model 
We use lower case to indicate 2D variables and upper case to 
indicate 3D variables. In 2D, the contact between two droplets 
has a contact length Uj, and in 3D, the contact has contact area 
Ajj. The force law for the two models are 



2D Model: fij = 72/) — , where r/y = — 



(2) 
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3D Model: Fi^ 



73D-^, where Rij 



RiRj 



(3) 



In the above equations, Y2D is a 2D Hne tension and y^o is a 
3D surface tension. For scenarios where droplets are asym- 
metrically deformed in 3D, the radius of curvature Rij in the 
3D model must be replaced by the mean curvature Rij. 

The 2D model would be straightforward to apply as we di- 
rectly measure lij, Vi, and Vj. To apply the 3D model, a rea- 
sonable assumption is that Aij is related to Uj and perhaps the 
droplet radii. The radii Ri and Rj are measurable as described 
in the previous subsection. 

Rather than choosing between the 2D and 3D models, we 
test generalizations of both models and let the data select what 
works best. As described above, one of our variables the force 
will depend on is Mij and we constrain all possible force laws 
so that f{Mij = 0) = 0. In general, we consider models of the 

form f^f^ {Alij, l/nj; a) for 2D and f^f^ (A/,-,-, I /Rij; a) for 
3D. a — ai,a2,... are the fitting parameters associated with 
a given functional form. To give an example, we could write 
fT^P — ai{Mij / Vij)^'^ with fitting parameters a\ and a2. In all, 
we test a total of 86 various 2D and 3D force laws of different 
functional forms that include exponentials, hertzians, power 
laws, and polynomials in Uj, l/^j, and I /Rij, and combina- 
tions of these forms. 

A.3 Optimization Problem 

To test the force laws, we establish constraints from the data, 
optimize each force law subject to the constraints, and then 
quantify how well the optimum force laws describe the data. 
To start with, we consider the constraints on forces in the x 
and y directions. 

In the _y-direction the sum of the forces on any given droplet 
is equal to the buoyant weight Wd- This is in practice hard 
to use directly, as Wd is small compared to the contact forces, 
and likely below limits set by noise. Therefore, rather than 
considering individual droplets, we note that droplets located 
at a given y must support the observed total buoyant weight 
Wobs of droplets below them, known simply from measuring 
the total area of droplets with centers below y. The way in 
which these droplets support this buoyant weight is through 
contact forces, and for an assumed force law fij{Alij^ 1 / rij; a) 
we can determine these contact forces by substituting our mea- 
sured values for Alij and rij (or Rij) into the function. If the 
assumed force model accurately predicts the forces, then the 
sum of these contact forces Y^^modj at a given y will equal 
Wobs- Here Y.^mod,y are the sum of the j-component of only 
those forces pointed in the downward directions. The reason 
we only consider the downward facing forces is because the 
collective buoyant weight is pushing upward, and to satisfy 
Newton's 3rd law, the balancing forces must be facing down- 
ward. We convert Wobs and F^od^y into 2D pressures (force per 



unit length) by writing Xobs = ^obs/^-> Kiod = "LJ^modj/^. us- 
ing the width of the chamber w. A is in essence the 2D hydro- 
static pressure at height y. Because there Js no static friction 
at the sidewalls, there is no Janssen effect^^. 

We define a goodness of comparison in the j-direction as 



^ [{^{y)obs -^{y)mod) / {^{y)obs)]^ , (4) 



where smaller values of Xy indicate a better match between 
the assumed force law and the actual forces. In the equation, 
y indexes various distances up the incline where )^{y)mod and 
^{y)obs are sampled, and the angle brackets are an average 
over _y. We normalize by {X{y)obs) to make Xy dimensionless, 
and since {X{y)obs) is independent of the assumed force law, 
it does not change the results. We sample A at intervals of 5ro 
up the incline. At each y sampled, A^^j is calculated using 
the contact lengths and droplet radii for all droplets found be- 
tween a distance y — Sr^ and j + 5ro up the incline, and Xobs is 
calculated using the position and radii of all droplets below a 
distance y up the incline. 

We next consider the forces in the x-direction. In contrast to 
the j-direction there are no external forces, so the sum of the 
forces on each droplet in the x-direction is zero. From this we 
construct the goodness of comparison 



ZUij /{\fi\) 



(5) 



where the F^^ij is the x component of the force at a contact 
between droplets / and j and (|7;|) is the average net contact 
force exerted on droplet /. In the equation, f^^ij are the forces 
predicted by the assumed force law. Due to measurement er- 
ror, the forces will not sum to zero, and the deviation from zero 
grows with (|7^|). We assume that the deviation will grow lin- 
early with (1 7; I) and to fairly weight the contributions of each 
droplet to xl^ we normalize the sum of the forces by (|7; |). 

Finally, we define a net goodness of comparison x^ = xlXy 
which indicates how well an assumed force law models the 
forces in both the x and y directions. Since we know the buoy- 
ant weight of our droplets in units of /iN, this allows us to find 
a force law in units of /iN. Later, we compare x^ between the 
different force laws to determine the best overall force law. 

A.4 Empirical Force Law: Monodisperse and Bidisperse 

We now apply our method to find an empirical force law. We 
start by determining the force law for same size droplets in 
contact using data taken on four different monodisperse sam- 
ples. The samples are prepared by placing droplets with 3D 
radius Rq into a sample chamber with gap thickness either 
100 ±4 jim or 180 ±4 /im, and once in the chamber, the 
droplets have a 2D radius of tq. The error assigned in the 
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Ro [Aim] 


ro [Aim] 


h [/im] 


poly (%) 


h [Mm] 


164 


183 


186 


5.5 


50.8 


143 


156 


180 


1.6 


48.8 


105 


128 


96 


1.9 


36.6 


84.1 


89 


96 


2.5 


30.3 



Table 3 Parameters characterizing the droplets in our 4 different 
monodisperse samples. Rq is the 3D radius of the droplets, tq is the 
2D radius, h is the gap thickness of the chamber, poly is the 
polydispersity of the sample, and /q is the length of contact for two 
droplets just in contact. Our uncertainties for the various measures 
are ± 0.2 fdm in Rq, ± 2 /im in tq, and ± 4 /im in h. The variability 
in Iq grows with droplet size and the measurement uncertainty can 
be expressed as ± 0.04 ro in Iq. 



gap thickness represents the unavoidable variations in the gap 
thickness measured at different points along the sample cham- 
ber. We note that the larger the sample chamber, the larger 
the variability of gap thickness we observe. As the chambers 
are inexpensive, sample chambers with more variability than 
±4 Aim are discarded. After the sample chambers are filled, 
they are sealed to prevent evaporation, and then placed on a 
microscope inclined at 28°. Droplets rise to the top and come 
to rest in mechanical equilibrium, at which point we acquire 
images of the sample. Various parameters characterizing each 
monodisperse sample are shown in Table [3] 

To determine the best force law we pick each possible func- 
tional form, optimize the parameters a for it, and then com- 
pare for the different functions. Several functional forms 
all have small values, and of these we choose one that is 
simple and plausible. For simplicity, measuring the 2D rij is 
simpler than measuring the 3D Rij. For plausibility, functions 
that treat Alfj and rij as ^kj/rij most closely resemble Eqn.[3] 

We judge the most reasonable function with low x^ ^ be 
Fij = aiAlij/rtj + a2{Alij/rij)^, quite similar to Eqn. [3] We 
also choose this function because we find that we can easily 
rescale this function using the 3D droplet radius Rq so that 
we have a universal force law for all four data sets (the four 
different droplet sizes). The rescaled force law is 



Fij = aiRoAlij/nj + a2{RoAlij/rij) 



(6) 



where ai = 4.25 AiN/mm and a2 4.12 jiWmrn^. 

The rescaled force law is shown for each data set in Fig. [13 
We see that all the force laws are nearly linear; the quadratic 
correction is about 10% for the largest forces. The rescaled 
force law shows that all the data collapse very well and only 
slightly deviate between each other at larger Alij/vij. For these 
larger values in AUj/ rij, the area fraction is close to 0.96 which 
is the upper limit where we can still confidently measure rij; 
for those close-packed droplets, only a few pixels occupy the 
water-oil interface and r,,- is hard to determine. We do not 
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Fig. 13 (Color online). A plot of the universal empirical force law 
for each data set. The solid lines are samples with a gap thickness of 
approximately 180 Aim and the dashed lines are samples with a gap 
thickness of approximately 110 Aim. The legend indicates Rq. 



claim that Eqn. [6] is the correct force law; in particular, while 
a\ has units of surface tension and is plausible for an oil-water 
surface tension, the physical meaning of is unclear. Rather, 
Eqn. |6] accurately provides the forces between our droplets, 
within the measurement limitations set by our data. Also, 
there may be other sources of error, for instance, Lacasse et 
al. has numerically shown that the force law has a slight 
sensitivity to the number of neighbors and the relative posi- 
tioning of the neighboring droplets. To examine if there are 
other potential sources of error, using Eqn. [6] we compared the 
deviations in the computed net force on each droplet to the de- 
viations we expect given our measurement errors, and find the 
two agree well. Thus, within the limitations of our measure- 
ment errors, we have resolved the forces as best as possible, 
confirming Eqn.|6]is adequate. 

To test how well one can determine a force law given fi- 
nite data and measure error, we additionally simulated inclined 
mechanically stable droplet packings of 1000 droplets with 
a known force law, then added noise to the data consistent 
with experimental noise. Applying our empirical method to 
the simulated data, we recover the known force law with 2% 
errors in the coefficients (noise equivalent to the experiments 
using the 5 x microscope objective) or 5% errors in the coeffi- 
cients (noise equivalent to the experiments with the 1.6x mi- 
croscope objective). This suggests it is possible that the ~ 4% 
variations between the force laws for different sized droplets 
seen in Fig. [13] are simply due to noise, and they may well 
have exactly the same force law. 

So far we have focused on force laws in monodisperse sam- 
ples, but we also need to measure forces between different- 
sized droplets in bidisperse samples. To obtain a force law 
between droplets of different sizes, we apply our method to 
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Table 4 The parameters characterizing the droplets in our 3 
bidisperse samples. Other parameters related to these droplets are 
shown in Table 1 The first column a = r^^^^^ jr^^^^^^^ is the size 
ratio. The two radii shown are the 3D droplet radii of the small and 
big droplets before placing them in the chamber. There are also 
three contact types: small-small {ss), small-big {sb), and big-big 
{bb), and therefore there are three /q values. Our uncertainties for 
the various measures are ± 0.2 Aim in /?o and ± 4 Aim in h. The 
variability in /q grows with droplet size and the measurement 
uncertainty can be expressed as ± 0.04 rg. 



with moving droplets, overlapping images of different fields 
of view would be difficult or impossible. This situation would 
require limiting the field of view to fewer droplets, if the 
same resolution of forces was desired. For any magnification, 
uniformity of lighting is essential so that the appearance of 
droplets is uniformly related to their true shape and size. As 
discussed in the previous subsection, our imperfect knowledge 
of the force law gives us a systematic uncertainty no worse 
than 5%. To determine the random uncertainty for particular 
forces, we take measured Vij and Mij values, add noise com- 
mensurate to our known uncertainty (discussed in Sec. |A.1| ), 
and recalculate the force to see the variation. The bidisperse 
data of Sec. |4] were taken with a 1.6x lens and have a ran- 
dom uncertainty of 16%. The monodisperse data used in the 
calibration procedure were taken with a 5x lens and have a 
random uncertainty of 8%. 



find an empirical force law using data taken on three differ- 
ent bidisperse samples. The bidisperse samples are prepared 
in the same manner as the monodisperse case. Table |4] sum- 
marizes the various parameters of our bidisperse systems; see 
also Table [T] 

For the case of a bidisperse sample with small and big 
droplets, there are 3 possible contact types to consider: small- 
small, small-big, and big-big. Our previous results give us 
small-small and big-big forces. We assume the unknown 
small-big force law obeys the same functional form as the 
monodisperse case (Eqn J6]), where a\ and a2 need to be deter- 
mined. Recall that Eqn.|6|contains a term Rq that rescales the 
force law and makes it universal. For the small-big contacts 
there are two different R^ values, one for each droplet size. To 
account for these two radii we substitute R^ with the arithmetic 
mean of the two radii (7?o) giving as our bidisperse empirical 
force law Fij = ai (Rq) Alt j/rij + a2{{Ro) Ali j / rt j)^ , where ai 
and a2 are unknown. To obtain ai and a2 for our bidisperse 
samples we minimize Xx^ ^nd find that ai and a2 are very 
close to that found for the monodisperse case and within the 
4% variation we expect from due to finite sampling and mea- 
surement error. Therefore, we have shown that to within 5% 
error we have a found a universal force law that works for any 
droplet size and is close to Princen's 2D with a 

small second order correction. 

A.5 Final comments on force law 

The uncertainties in determining forces are related to the mag- 
nification. The higher the magnification, the better we can 
measure the contact length / and the mean curvature r. For- 
tunately, given that we study static samples, this means we 
can take overlapping images at high magnification to reduce 
our uncertainties, as described in Sec. |2] In an experiment 



An additional experimental complication is that droplets at 
rest occasionally feel a static force from the top and bottom 
plates. This is likely due to contact line pinning on impurities 
or microscopic scratches on the glass. To minimize this, we 
pre-clean each slide with methanol which we gently blow off 
the slide. Harsher cleaning methods do not significantly re- 
duce the droplet pinning. The magnitude of these forces can 
be estimated by examining a dilute concentration of droplets 
in a horizontal slide, and then slowly tilting the slide to see 
when the droplets begin to move due to gravity. For the sam- 
ples discussed in this work, they begin to move at tilt angles of 
about 4.5° or sooner. We discard any sample chambers with 
pinning stronger than this. Given that our experiments are con- 
ducted at a tilt angle of 28°, the buoyant weight of a droplet 
is sin(28°)/ sin(4.5°) = 6 times larger than any pinning force. 
Taking the analysis a little further, since the friction force on 
any droplet can range from zero to the maximum, a more 
appropriate estimate for the buoyant weight of a droplet is 
sin(28°)/(sin(4.5°)/2) = 12 times the average pinning force. 
These estimates show that the inter-droplet forces seen in the 
jammed emulsions ((j) > (\)c) are on the order of a hundred 
times the pinning forces. We believe that the situation in our 
calibration experiments are even more favorable. After com- 
paction, the pinning forces should be in random directions, 
as an analogous granular experiment observed that particles 
move in random directions during compaction Therefore, 
a vector average gives a pinning force on each droplet very 
close to zero. Since our empirical method relies on the aver- 
age vector forces on a droplet, the pinning forces can safely be 
neglected in the results of Sec.|4] Note that if our experiment 
was scaled down in size (smaller droplets, thinner plate gap) 
the pinning forces become more significant compared to the 
droplet weight and can dominate the results. 
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